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•^ ■ Abstract 

f"^ I Attention has been paid to the similarity and duality between the Gregory-Laflamme instability 

of black strings and the Raylcigh-Plateau instability of extended fluids. In this paper, wc derive a 
set of simple (l+l)-dimcnsional equations from the Navier-Stokes equations describing thin flows 
f^ , of (non-relativistic and incompressible) viscous fluids. This formulation, a generalization of the 

O I theory of drop formation by Eggcrs and his collaborators, would make it possible to examine 



^ 



X 



the final fate of Rayleigh-Plateau instability, its dimensional dependence, and possible self-similar 
behaviors before and after the drop formation, in the context of fluid/gravity correspondence. 
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1 Introduction 

Black holes in higher dimensions have been studied from various points of view. In particular, they 
play practical and important roles in the contexts of large extra-dimension scenario and gauge/gravity 
correspondence. It has become known that they have rather richer structures than reported. Some 
explicit examples of instability and non-uniqueness in a solution space have been presented. In such 
a case, the phase transitions among various phases of black holes are expected (see [T] for a review). 

One important class of higher-dimensional black holes is that of Kaluza-Klein or 'caged' black 
holes, i.e., black holes in a spacetime one of whose spatial dimensions is compactified on a circle, S^. 
Black objects known in this system include uniform black strings, non-uniform black strings, and black 
holes. Here, the black strings wrap the S^, whereas the black holes are localized in the S^-direction, 
and 'uniform' stands for the translational invariance in the S^-direction. See [21 [3] for reviews. The 
uniform black strings were shown to be unstable by Gregory and Laflamme [3] with respect to linear 
perturbations breaking the translational invariance. Furthermore, the uniqueness is explicitly broken 
and the transitions among the three phases are expected. It is noted that the transition from a black 
string to a black hole via the pinch-off of horizon is accompanied by a topology change. Horowitz and 
Maeda [5], however, argued that the Gregory-Laflamme instability does not result in the pinch-off of 
horizon in a finite affine time. After this, a numerical simulation was performed by Choptuik et a/. [B], 
and many studies on the sequence of static solutions were done to figure out the phase diagrams HI 

The gauge/gravity or AdS/CFT (anti-de Sitter/conformal field theory) correspondence has pro- 
vided the powerful tools to investigate strongly coupled quantum field theories by examining classical 
gravity. The behaviors of the field theories at finite temperature are encoded in the thermodynamics 
of black holes residing in higher dimensions. Thus, from the viewpoint of black hole physics, it is 
natural to expect that we could learn about the black holes from the field theory side. Regarding this 
point, there was important progress called the fiuid/gravity correspondence. It was shown that the 
Einstein equations with a negative cosmological constant are dual to the Navier-Stokes equations for 
a certain kind of conformal fiuid residing on the AdS boundary [8] (and see [9] for a review). This 
seems to open a new way to examine problems in black hole physics by working on the fluid mechanics 
which is a collective description of field theory at a long wavelength limit. 

The similarities between the Gregory-Lafiamme instability of black strings and the Rayleigh- 
Plateau instability of fluids supported by surface tension or Newtonian gravity were pointed out 
by Cardoso and his collaborators \10\ lll|. and then the similarities were found to persist indeed up 
to non-linear regimes [12], at least as to the bifurcation structures in a phase diagram. One possi- 
bility of the physical origin of similarity is that the Gregory-Laflamme is holographically dual to the 
Rayleigh-Plateau instability of fluid lumps supported by surface tension, which is suggested by the 
conjecture that a lump of 'gluon' plasma in large- A/" gauge theories is dual to a black hole localized in 
an asymptotically AdS space with a Scherk-Schwarz compactiflcation [131 EU • It is natural to expect 



^Recently, Lehner and Pretorius reported the results of a new numerical simulation of a 5-dimensional black string [Tj- 
They argue that the horizon indeed breaks up in a finite asymptotic time, accompanied by a self-similar cascade of the 
instability. 



that in the theories aUowing a first-order confine-deconfine phase transition the deconfined phase ap- 
pears as locahzed bahs separated from the confined phase by a certain kind of boundaries. The point 
of the above conjecture is that a gravitational dual of such a plasma ball is a finite-energy black hole 
(rather than an infinite-energy black brane) in the bulk (see Introduction of [15j for a brief but nice 
review). The phase diagrams of plasma tubes and plasma balls obtained in |16j and [T7] with this 
approach mimic well those of the Kaluza-Klein black holes. These results suggest that we will be able 
to examine the dynamics of Kaluza-Klein black holes in non-linear regimes by working on the fluid 
mechanics, though the duality originally proposed in [13l HI] is still a conjecturecl 

The fluid equations involved in the fluid/gravity correspondence are the relativistic Navier-Stokes 
equations with speciflc values of transport coefficients (viscosities, thermal conductivity etc.) and an 
equation of state, that are much easier to treat than the dual gravitational equations. When the 
non-linear dynamics of a fluid is concerned, however, the fluid system is still difficult to treat as it 
is. Therefore, it is of practical importance to import and/or develop techniques in fluid mechanics to 
investigate the problems in black hole physics via the correspondence. In this paper, we import the 
scheme developed by Eggers and Dupont [20] to describe the drop formation of free-surface flows (see 
|2HI22j for comprehensive reviews), while generalizing it to that in arbitrary dimensions. The essence 
of this scheme is to reduce the (2 + l)-dimensional partial differential equations and their boundary 
conditions describing axially symmetric flows to a set of (1 + l)-dimensional equations by a radial- 
expansion method which is valid when the characteristic length scale in the direction perpendicular to 
the axis is smaller than that in the axial directiono It is noted that the approximate equations become 
exact as the pinch-off of a fluid neck is approached. In addition, we will see that with an appropriate 
prescription the approximate equations can describe relatively thick flows such as those at the onset 
of instability. By the vanishing of one independent variable, the problem is significantly simplified. 
Furthermore, in this approach the boundary conditions describing the balance among the surface 
tension, viscosity, and pressure jump at the fluid surface are contained in the set of reduced equations. 
In other words, the boundary conditions 'disappear' from the system. While in this paper we will 
focus on a non-relativistic and incompressible limit for simplicity, the generalizations to relativistic 
and compressible cases would be straightforward (but see, e.g., [231 IMl [25] for non-relativistic and 
incompressible fluids obtained via the fluid/gravity correspondence or the membrane paradigm). 

In Sec. [21 we introduce the Navier-Stokes equations and their boundary conditions, and the reduc- 
tion to the axially symmetric flows is made. In Sec. [31 the Navier-Stokes equations and their boundary 
conditions are combined to derive the set of (1 + l)-dimensional equations, and their properties are 
investigated. In particular, the growth rate of the Rayleigh-Plateau instability of viscous flows are 
derived. Sec. [His devoted to discussions. 



■^Another powerful and efficient way to study tfie Gregory-Laflamme instabiUty is to solve the Einstein equations for- 
mulated in terms of an effective fluid that lives on a dynamical worldvolume, which is called the 'blackfold' approach |18) . 
The long- wavelength component of Gregory-Laflamme instability was obtained as a sound-mode instability [181 119| . It 
is noted that this approach does not need a boundary of fluid, and so the instability is not Rayleigh-Plateau's one. 

^Needless to say, A*' -I- 1 means A'' spatial dimensions plus time, ft is noted, however, that descriptions by (A^ + 1)- 
dimensional partial differential equations is said A'-dimensional (as the title of this paper) in non-relativistic physics. 



2 Viscous flows with boundaries 

We review the Navier-Stokes equations and their boundary conditions (see, e.g., [26]), while gener- 
ahzing them into those in arbitrary dimensions and curvihnear coordinates. See |141 \T7\ [T5] for a 
relativistic formulation and especially [12] for a viscous case. 

2.1 Navier-Stokes equations 

We consider a fluid in a d-dimensional flat spacetime {d > 4) with a time coordinate t and general 
curvilinear spatial coordinates x^ , (/, J = 1, 2, . . . , d — 1). The equation of continuity may be written 
as 

dtp + Viipv^) = 0, (2.1) 

where p is the mass density, v the fluid velocity, V/ a covariant derivative with respect to spatial 
metric gjj{x ). The equations of motion for the fluid are obtained from the momentum conservation, 

dt{pv')+VjU'-^ = 0, (2.2) 

where II^'^ is a momentum flux density tensor. The momentum flux density tensor in general may be 
written as 

II/j = pvivj - aij , (2.3) 

ajj being a stress tensor. We consider the viscous fluid whose stress tensor is given by 

o"/j = -pgij + V ( ^ivj + Vjvi - -j—^aijO j + QgijO , (2-4) 

where p is the pressure; rj and ( are the shear and bulk viscosities, respectively, assumed to be positive 
constants; 6 := Vjv^ is the expansion. Plugging the above expression of momentum flux density tensor 
into the conservation law (j2.2p . we obtain the Navier-Stokes equations with the help of equation of 
continuity (|2.ip . 

pDtVi = -Vjp + 7]Avj +U + ^^^j ^i^ ' (2-5) 

where the so-called substantial time derivative and Laplacian are defined by Dt = dt + f ^V/ and 
A = V'^Vj, respectively. 

Now, let us move to the boundary conditions imposed on the above equations when one considers 
two immiscible fluids (call them fluid 1 and fluid 2) divided by a surface. The position of the surface 
may be identified by the vanishing of a scalar function, f{t,x^) = 0. Denoting the unit normal of the 
surface by rij = V///|V/| (pointing from fluid 1 to fluid 2) and surface tension by a{t,x^), the stress 
balance equations, called the Young-Laplace formula, are given by 



.(2) _ ^(1)^ ^j ^ ^^^^ ^ P/Vja 
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(2.6) 
/=o 



where k is [{d — 2)-times] the mean curvature of the surface and Pjj is the projection tensor, 

d-2 

K = Y, K^ = V/n^ , PiJ = giJ - ninj . (2.7) 

Here, Ri is the curvature radius (assumed to be positive if it is drawn in the side of fluid 1) in a 
principal direction on the surface. 

Here, we assume a few properties of the fluid and surface for simplicity. In the rest of this paper, 
we will focus on an incompressible fluid {Dtp = 0), for which equation of continuity (j2.ip and Navier- 
Stokes equations ()2.5p simplify, 

Viv^ = , (2.8) 

Dtvi = -^ + uAvi , (2.9) 

^ = r]/p being the kinetic viscosity. Second, we assume that surface tension a is constant and neglect 
the effects of fluid 2 (i.e., ajj = 0), in which case the normal and tangential components of stress- 
balance equation (j2.6p are given by (omitting the label of fluid 1) 

—nja nj = an L^„ , (2-10) 

P/ajKU^ = l/=o ' (2-11) 
where the stress tensor is reduced to 

(TIJ = -pgij + ri{Vivj + Vjvi) . (2.12) 

Finally, we present a kinetic boundary condition called the free-surface condition, stating that the 
surface moves with the fluid, 

Dtf = l^^o . (2.13) 

Hereafter, we will consider Eqs. (J2.8p - (j2.13j) . 

Under the above assumptions {i.e., the incompressibility, and the constancy of surface tension and 
viscosities), free parameters of the fluid are the density p, kinetic viscosity v, and surface tension a, 
which are dimensionally independent each other. With these parameters, one can define typical length 
and time scales ac| 

£.:=^, t.:=f^. (2.14) 

a a^ 

It is noted that one can regard, for example, {p,iu,tu) as a set of three independent fluid constants 
rather than {p, u, a) by reversing (I2.14p , 

9^, . = 'f. (2.15) 

''1/ ^v 




Figure 1: A schematic figure of the axially symmetric flow in (n + 2)-dimensional space. 



2.2 Axially symmetric flows 

Since we are concerned with the dynamics of axially symmetric flows, it is convenient to introduce a 
cylindrical coordinate system with which the line element of {d — l)-dimensional flat space is written 
as 



gij{x^)dx^dx-^ = dz^ + dr^ + rMf^^ = Sabdx^dx" + r^7ij((^*)d<^M0J , 



(2.16) 



where dQ^ = 7jj((/))d(/)*d(/>-^ {i,j = 1,2, ... ,n := (i — 3) is the line element of a unit n-sphere and 
x" = {z,r). Now, let us assume that the fluid has a velocity only in the directions of z and r, and 
that the fluid surface respects the axial symmetry. That is, the velocity field v^{t,x^) and the scalar 
function f{t,x^) are written as 



v" = v''{t,z,r) , v' = 0, f = r-h{t,z), 

where h{t, z) is called a height function, giving the local radius of flow (see Fig. [T|). 
Then, equation of continuity (|2.8p reduces to 

dzVz + (dr + -)vr = . 



(2.17) 



z and r components of Navier-Stokes equation (j2.9p are given bjo 



OtVz + {Vzdz + Vrdr)Vz 



+ V [dl + dl + -d^ Vz , 

dtVr + {Vzdz + Vrdr)Vr = — + U ( d"^ + d'^ + -dr ^ 

n \ J. j.^ 



dzP 

P 
drP 



Vr 



(2.18) 

(2.21) 
(2.22) 



'^Iv and iy may be called viscosity length and viscosity time, respectively. Their order of magnitude considerably 
depends on the kind of fluids (and interfaces), e.g., (.-^ ~ 10~* cm for mercury and tv ~ 10^ cm for golden syrup [22] . 



^Non-zero components of Levi-Civita connection F^k are 



pa 
J- ij 



-rSr^ij , Vl^ 



515. 






(7) pi 



- jk 



(2.19) 



where '^'F*^ is the connection with respect to the metric of unit n-sphere. With using these connections, the vector 
derivatives in the Navier-Stokes equation are calculated. For example, 



ViVj = daVb5'i5j + rVr'yij5}5j , AVa = {d^ + dr)Va H drVa J^allr 



(2.20) 



Normal and tangential stress-balance equations ()2.10p and ()2.1ip reduce to 



2v 



P 



P 1 + {d,hf 



{dzhfdzVz - {dzh){dzVr + drV^) + Sr^r 



1 + {d,hY 
where the mean curvature is given by 



2{dzh){dzVz - drVr) - [l - [dzhf] {dzVr + dr 



a 
—I 

P 




r=h 



r=h 



dlh 



+ 



n 



[1 + (a,/l)2]3/2 r,/TT{dj^ 



(2.23) 
(2.24) 

(2.25) 



It is noted that the z and r components of tangential stress-balance equation (j2.1ip lead to the same 
result, Eq. (|2.24p . due to the isotropy of surface tension. Finally, kinetic boundary condition (I2.13P 
reduces to 



dth — Vzdzh + Vr = 



3 One-dimensional approximation 



r=h 



(2.26) 



In this section, we derive a set of (1 + l)-diixLensional equations from the (2 + l)-dimensional equations 
obtained in Sec. 12.21 while generalizing the formulation in j20j (details are given in j21j). 

3.1 Radial expansion 

Let us consider the situation in which a characteristic length scale in the r-direction (typically, the 
thickness of the fluid neck) , denoted by ir , is smaller enough than a characteristic length scale in the 
axial direction, denoted by iz- We expect in this situation that the dynamics is dominated by the 
motion in the z-direction. Assuming this situation, we define a small dimensionless parameter e and 
a typical time scale r determined by the motion in the z-direction, 

e~-^, T~-^. (3.1) 

As will be seen below, we can show that the length and time scales, iz and r, in such a situation have 
the following order of magnitudes measured in the fluid constants (j2.14p , 



Iz ^ eiy , T ^ e ty , 



(3.2) 



if one demands the balance among the contributions to the dynamics from the inertia, pressure gra- 
dient, viscosity, and surface tension. 

First, we demand the balance among the three contributions to Navier-Stokes equation ()2.2ip from 
the inertia, pressure gradient, and viscosity: 



OtVz ~ ~ J^O^Vz 



(3.3) 



Replacing the derivatives in Eq. ()3.3p by corresponding scales, we have two relations for l^ and r, 

£2 p e 



~ u 



r^ p 



(3.4) 



One can estimate the order of pressure by the equation of normal-stress balance (|2.23p , namely p ~ an. 
In addition, the mean curvature of a sufficiently flat surface is dominated by the curvature in the 
spherical part [i.e., the second term on the right-hand side of Eq. (|2.25|) ] to yield k ~ l/^. ~ l/{e^z)- 
Thus, we have 



a 

If. 



(3.5) 



Eliminating the pressure in (j3.4p by (j3.5p . one can immediately solve (j3.4p for Iz and r, to obtain 
Eq. daj). 

From the above argument, we can estimate the order of magnitudes of all quantities in our basic 
equations in Sec. 12. 2i Namely, the magnitude of each quantity can be expressed in terms of dimen- 
sionless parameter e and fluid constants {p,iu,t,y). With this in mind, we introduce the following 
dimensionless variables with hats. 



t = e tyi , z = elyZ , r = e i^f , h = e £yh 



iu 



p^l 



— Vz , Vr = —Vr , P- 2^2 



P , K 



eH,. 



(3.6) 



Plugging Eqs. (j2.15p and (j3.6p into the basic equations, all equations are made dimensionless. 
Equation of continuity (|2.18p does not change its form. 



dzVz + (dr + -]Vr = ^ ■ 



Navier-Stokes equations (I2.21J) and (I2.22P lead to 

dtVz + {Vzdz + Vrdr)Vz = 
dtVr + {Vzdz + Vrdr)Vr = 



dzP + 

1 



5^ + 4 (dl + Vdr 



2 y^r 



r 



drP + 



gZ Y J, J. 



Normal and tangential stress-balance equations (j2.23p and (j2.24p become 



i + €^{dzhY 



e^{dzhfdzVz - (dzh) (e^dzVr + drVz) + drV, 



f=h 



1 



1 + e^dzh)^ L 
where the mean curvature is given by 



2e\dzh){dzVz - 4t)r) - 1 - e\dzhf {e'^dzVr + drVz) 



(3.7) 

(3.8) 
(3.9) 

(3.10) 
(3.11) 



e^dlh 



3/2 



+ 



n 



l + e'^{dzhf fJl + e'^{dzhf 



(3.12) 



Finally, kinetic boundary condition ()2.13p does not change its form 

— dth — Vzdzh + Vr = , . (3.13) 

f=h 

Now, we are ready to derive an effective (1 + l)-dimensional equations by approximating the 
dependence of fluid quantities on f by Taylor-series expansion. First, the velocity in the z-direction 
and the pressure are expanded as 

Vz{i,z,f) = vo{i,z)+V2{i,z){eff + 0{{er)^) , 
p{i,z,f) = po{i,z)+p2{i,z){erf + 0{{ef)^) . (3.14) 

Then, the velocity in the r-direction is determined order-by-order from the equation of continuity 
(1231), 

evr{i,z,f) = ^—dM^r) ^—dMerf + 0{ierf) . (3.15) 

n + 1 n + 6 



The leading order of Navier-Stokes equation in the r-direction (j3.9p gives 

4^0 + vo^t'o = -dzPo + 2(n + l)i)2 + dzVo ■ (3.16) 

The leading orders of normal and tangential stress-balance equations, (I3.10p and (j3.1ip . give 

2 

Po H —rdzVo = K , (3.17) 

n + 1 

2(n + 2)dzhdzVo + hd'^vo - 2(n + l)hv2 = , (3.18) 

where the leading order of dimensionless mean curvature is given by 

k = ^ . (3.19) 

h 

Eliminating pQ and V2 from Eq. (j3.16p with using Eqs. (j3.17p and (j3.18p . we obtain an equation of 
motion for vq, 

dtvo + vodzVQ = -dzK + ' / — - ,, ■ (3.20) 

Finally, from the leading order of kinetic boundary condition (j3.13p . we obtain an equation of motion 
for h, 

dth + ijodzh + hdzvo = . (3.21) 

n + 1 

We have obtained a closed set of (l-l-l)-dimensional equations, (|3.20p and ()3.2ip . that describes the 
flows of viscous fluid in general dimensions. Apart from the simplification that independent variable 
r is eliminated from the system, the leading-order effect of the stress-balance equations at the surface 
is encoded in the first term of the right-hand side of Eq. (|3.20p and the stress-balance equations 
themselves apparently disappear from the system. 



3.2 Conservation laws 

Let us write Eqs. ()3.20p and (j3.2ip in dimensionful form. Writing the corresponding dimensionful 
quantity without a hat and putting v{t,z) := VQ{t,z), we obtain a set of equations for v{t,z) and 
h{t,z), 

v + vv' = + — ^ , ^/ , 3.22) 

p n+1 /I'^+i ' ^ ^ 

h + vh' + hv' = , (3.23) 

n + 1 ^ ^ 

where the dot and prime denote the derivatives with respect to t and z, respectively. From Eq. ()3.19p . 

the mean curvature is given by 

K = ^ . (3.24) 

h 

Although the greatest significance of the system (I3.22P and (13.23P lies in the asymptotic limit of 

e — 7- 0, it has been successfully used for simulations away from the breakup (e.g., in [20]). The key 

to this success lies in a modification of Eq. (j3.24p to the complete expression of mean curvature [see 

Eq. ^2B], 

h" n , , 

+ , • (3.25) 



(1 + /i/2)3/2 hVlTh!^ 

The right-hand side of Eq. ()3.25p includes an infinite sequence of terms of the radial expansion. The 
above prescription has not been justified completely so far and is 'phenomenological' at this point. It 
is highly desirable to develop a more general and consistent method to 'resum' all the relevant terms 
of the expansion. Readers interested in this point are directed to Sec. V.B of [21], where additional 
rationale for this prescription is discussed. Anyway, we will investigate the dynamical properties of the 
system described by Eqs. (j3.22p and (I3.23P with the complete expression of mean curvature, Eq. (I3.25p . 
We will see that the above replacement indeed broadens the applicability of Eqs. (]3.22p and (|3.23p . 
In order to see the structure of Eqs. (j3.22p and (|3.23p . it is instructive to rewrite them as follow, 

dt{h^+\) + d,{h^+^v^) = -— /i"+i + ^i!i±^i.(/i"+i«')' , (3.26) 

p n+1 

a,(/j"+i) + a,(/i"+ir;) = 0. (3.27) 

The form of Eq. (|3.26p is best motivated by considering the force balance on a slice of fluid which 
carries linear momentum per unit length h'^'^^v, disregarding the common factor pO„, where O^ := 
27r*-"'"'"^^'^/r[(n + l)/2] is the area of unit n-sphere. The left-hand side is the total time derivative of 
the momentum. The first term on the right-hand side comes from the capillary forcing on the slice, 
the second term is the viscosity forcing. The meaning of Eq. (|3.27p is also clear: it is written as a 
conservation law for the mass per unit length h'^'^^, disregarding the common factor pil„ again. 

Now, we express the mass conservation and energy dissipation due to the viscosity in integral form. 
We write the mass in a fixed interval z £ [z-, 2;+] as 

M{t) = p^n / hJ'^^dz , (3.28) 



Then, one can show that from Eq. ()3.23p this mass satisfies a conservation law, 

dtM = -pQn (h'^+^v) '^ . (3.29) 

The kinetic energy and the potential energy stored in the surface may be given by 

^kin(i) = -P^n / h'^+^v'^dz , ^surf(i) = (n + 1)J^„Q / h^'^/l + h'^dz . (3.30) 

Combining Eqs. (I3.22p . (13.231) . and (I3.25p . one can show 

dt{E^^ + i^surf) = - ^^''+^f- p^ r h-^\'^dz 



^ /I ,r,4-i ^ 2(n + 2) ,„,, , ,„,i , , /i"/i'/i 

-^n -ph^'^^v^ - ^ > i//t"+^W + ah^-^^vK - (n + 1)q- 



2^^ n + l ^ ' ' VI + /i'2 



(3.31) 



This means that apart from the driving through the boundary terms, the total energy of the system 
can only decrease due to the viscosity, and the fluid eventually has to reach some static equilibrium 
shape corresponding to a minimum of the surface energy. It is noted that Eq. (j3.25p is essential to 
derive Eq. (I3.31[i^l and to have the above physically reasonable picture on the final equilibrium, which 
will be realized far from the asymptotic limit of e — )• 0. 

3.3 Instability of viscous flows 

The dispersion relation or growth rate of Rayleigh-Plateau instability of viscous flows can be derived 
analytically with using the above formulation. 

Let us consider linear perturbations of a static cylindrical tube whose radius is rg, 

h{t, z) = ro [1 + Se'^^ cos{kz)] , v{t, z) = ^Vbe""* sin(A;z) , (3.33) 

where 6 and Vq are constants {\5\ ^ 1). Plugging these expansions into Eqs. (j3.22p and (|3.23p . we 
obtain algebraic equations for w and Vq at 0{5), which are solved to yield 



u{k) = ^uo (Jj^^^ikro? [n - {kro)'] + ^kr^f - \^{kroA , (3.34) 

n + l \V {n + 2Y ro y tq / 



The growth rate in the low-viscosity {iu/fo ^ 1) and high- viscosity {iu/i^o ^ 1) limits are given by 

/ \ 1/2 



V"+l 



(A:ro)\/n-(A:ro)2 and ^- 2(^ + 2) (/) h - (^^o)^] , (3.36) 



6i 



We have used an identity, 
which holds for Eq. (IT^ but not for Eq. ^^. 
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(3.32) 



respectively. Both the above hmits for n = 1 coincide with the classic result of Chandrasekhar j27j if 
an expansion to lowest order in /erg is made. 

One sees from Eq. (j3.34p that the wavenumber of the onset of instability [kc > 0, uj{kc) = 0] is 
given by kc = y/n/rQ , irrespective of the viscosity, which generalizes the result of the inviscid case |10| 
and is similar to the dimensional dependence of Gregory-Laflamme instability [281 [29] . One the other 
hand, the wavenumber of the most unstable or fastest growing mode fcmax is a decreasing function of 
viscosity, 

{kraa^rof = — r 1 • (3.37) 



l + (n + 2)V(4Ao)/(n + l) 

Thus, in the high-viscosity limit a mode of infinitely long wavelength becomes the most unstable one. 
The dependence of the dispersion relation on the dimensions and viscosity are depicted in Fig. [2j One 
can see that the growth rate of the most unstable mode, a;(/cmax), is also a decreasing function of 
viscosity, which represents the slowing down of motion by the viscosity. For reference, a derivation of 
the growth rate for the inviscid case not relying on our formulation in Sec. [3] is given in appendix lAl 
It is interesting to see a long-wavelength limit with the viscosity kept arbitrarylj Up to the 
quadratic order of Zero (^ 1), growth rate (j3.34p becomes 



n + 1 



1 _ /^ + ^ /5 (kr ) 
^n(n + l)V ro ^ ''^' 



(3.38) 



Thus, one can see that the viscosity is responsible for the quadratic term, which is consistent with the 
expression for the low- viscosity limit in Eq. (j3.36p . It is noted that expression ()3.38p does not change 
if one uses the 'incomplete' expression for the mean curvature, Eq. ()3.24p . 

4 Discussion 

We have obtained the set of (1 -|- l)-dimensional equations, p.22p and ()3.23p . for the velocity in the 
axial direction v{t,z) and height function h{t,z), describing the dynamics of axially symmetric thin 
flows of non-relativistic and incompressible fluids, where the leading-order effects of the homogeneity 
in the radial direction are taken into account by series expansion p.l4p . In spite of the significant 
simplifications {i.e., the reduction of the numbers of independent variables from three to two, and 
the 'disappearance' of the boundary conditions resulting in the reduction of the numbers of equations 
from five to two), the set of equations have nice properties: it exhibits the mass conservation ()3.29p 
and the dissipation of total energy by the viscosity (|3.3ip in a plausible manner, provided the complete 
expression for the mean curvature (j3.25p is used. The growth rate of Rayleigh-Plateau instability (j3.34p 
was derived, producing the expected/reasonable dependence on the spacetime dimensions d = n + 3 
and shear viscosity r] = pv. 

This work has many directions of application and generalization, a few of which are 



^The author thanks R. Emparan for his suggestion to see this hmit. A relevant discussion on fc^-term can be found 
in Sec. 3 of Camps et al. [THj, though their approach to the Gregory-Laflamme instability itself is different from ours. 
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Figure 2: (a) The dimensional dependence of the growth rate of Rayleigh-Plateau instability for Ivjri^ = 1. 
d = 4, 5,6,7 from the bottom to the top. (b) The viscosity dependence for d = 5. iy/vQ = 0,0.1, 1, 10 
from the top to the bottom. 



Equations (|3.22p and (|3.23p with n = 1 has been used to describe successfully the drop formation 
|20j . In particular, it has been shown that there is a self-similar solution (precisely speaking, 
one is before the breakup and another is after the breakup) that plays the role of an attractor 
in the pinching region |301 [3T] . The analysis on higher-dimensional self-similar solutions will be 
reported elsewhere [32] rl 

One of interesting aspects of both the Gregory-Laflamme and Rayleigh-Plateau instabilities is 
their dimensional dependence of the phase structures [331 [M| and dynamics [35] that appears 
in non-linear regimes. Above a critical dimension, it is expected that there are non-uniform 
equilibrium configurations of black string and fluid flows that serve as the end point of the 
instabilities. Thus, it would be interesting to examine the dimensional dependence of the end 
point of Rayleigh-Plateau instability with Eqs. (j3.22p and (j3.23p rl 

It would be possible to generalize the formulation in this paper to the relativistic and compressible 
cases such as one in ]\?>\ [T^ since we would have no additional degrees of freedom essentially 
once an equation of state is provided by the fluid/gravity correspondence. 



*Such a self-similar solution describing the breakup takes the form of w ~ (io — i)~^ 1^(C) s^nd ft ~ (io — t)H{^) for 
d — 4. Here, {to,zo) is the time and location of the breakup, and H and V are regular functions of similarity variable 
^ = (2 — zo)/{to — i)^". It would be interesting to see whether the same behavior is found in higher dimensions. See [7] 
for the observation on the similar behavior near the breakup of black-string apparent horizon. 

^Due to the simplicity of this set of equations, rather 'economical' simulations would be possible. The time evolution 
of the instability can be traced, e.g., with the NDSolve command in Mathematica installed on a laptop computer, at least 
for simple {e.g., periodic) boundary conditions and away from singularities. 
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A Instability of inviscid flows 

Let us consider linear perturbations of an inviscid flow {v = 0) that has a uniform velocity in the r- 
direction. For this purpose, it is convenient to work in a boosted frame in which the uniform velocity 
vanishes. We assume that the perturbation results from a sinusoidal disturbance of the height function 
given by 

h{t, z) = ro[l + e{t) cos{kz)] , (A.l) 

where |e(i)| <^ 1. This disturbance of the height function leads to those of the pressure and velocity 
field, where the disturbed pressure would take form oi p = po + Sp{t, z, r). Here, po is a constant de- 
termined by the normal-stress balance as po = na/rQ. For incompressible fluids the pressure operated 
by Laplacian vanishes in general, Ap = 0, and therefore we have 

A6p = . (A.2) 

Writing the perturbation of pressure as 6p = 6p(t)Fk{r) cos{kz), Eq. ()A.2p leads to 

-r^ + --r^-k^Fk = o. A.3 

dr^ r ar 
This equation is solved by the modified Bessel function of the first kind {e.g., |36]). 

where we have assumed the regularity at the axis, r = 0. The combination of ()A.4p and the per- 
turbation of normal-stress balance ()2.23p (i.e. 6p\r=ro = (xSk) determines the time dependence of the 
pressure perturbation, 

I{n-i)/2{kro) 

The perturbation of velocity in the r-direction, 5vr{t, z, r), at the boundary is immediately obtained 
from kinetic boundary condition (12.260 . 

5vr\r=ro = dtSh = rQe{t) COS (kz) . (A-6) 

The Navier-Stokes equation in the r-direction, Eq. (I2.22p . relates the pressure- and velocity-perturbations, 

drSp 



dt6vr 



r=ro p 
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(A.7) 

r=ro 




Figure 3: The growth rate of Rayleigh-Plateau instability for an inviscid flow, in d = 4,5,6,7 from the 
bottom to the top. 



Plugging equations ()A.4p . ()A.5p . and ()A.6P into equation ()A.7p . we obtain an equation for e(t), 



e{t) = e{t) 



ar'^^-'^'^n-ikr.f]^ \I^n-i)/2{kr) 



pI{n-l)/2{krQ) 



r("-i)/2 



(A. 



r=ro 



After some calculations with assuming e{t) oc e"^*, we obtain a dispersion relation of the linear pertur- 
batiorE§ 



UJ = ijJl 



.kro[n - {krof]I(n+i)/2{krQ) 



'(n-l)/2 



{kro) 



(A.IO) 



which reproduces the result in \\\\ I17j. The dispersion relation for the inviscid case, Eq. (lA.lOp . in 
several dimensions are depicted in Fig. [3l 
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